hep-th/0506138 



A practical model for cosmic (p,q) superstrings. 



P. M. Saffin§ 

Department of Physics and Astronomy, University of Sussex 
Falmer, Brighton BN1 9QJ, UK 



Abstract 

We propose a model with a U(1)axU(1)b gauge symmetry which contains topological strings 
carrying magnetic flux under each of the U(l)s. By calculating the tension of the first few 
low energy strings we show that bound states, containing flux under each of the U(l)s, are 
stable against decay to their constituents of lower winding number. As the model contains only 
Abelian gauge symmetries this model provides a pragmatic solution to the numerical modelling 
of cosmic-superstrings. Whilst these defects do not satisfy a BPS bound we argue that they 
bear sufficient similarities to warrant such a study. 
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1 Introduction 



Cosmic strings have come in and out of fashion a number times, and whilst it is clear that they do 
not play a lead role in structure formation they may nevertheless have observational consequences 
El- Indeed, it is hoped that advanced LIGO and LISA will be able to detect the characteristic 
radiation of a string network [3J 03 10] • The context of the current revival of cosmic strings dates 
back to work of Witten in perturbative string theory [7] where it was concluded that superstrings 
could not act as cosmic strings. Since then it has become clear that there are more possibilities 
for the compact dimensions of string theory, which has led to the issue of cosmic superstrings 
being reconsidered [HJ 03 HUj ; for a nice overview see . The particular networks that we are 
interested in are associated to the low energy dynamics of type IIB string theory which, in ten 
dimensions, contains two different species of string. First of all there is the fundamental, F-string 
which carries charge under the Neveu Schwartz-Neveu Schwartz two-form potential, secondly there 
is the Dirichlet, D-string carrying charge under the Ramond-Ramond two-form potential. The types 
of string that motivate the current work are (p, q) strings which are a supersymmetric bound state 
of F and D strings. These come about because parallel F and D strings are not supersymmetric and 
can reduce their energy by the F string breaking on the D string with the endpoints moving off to 
infinity |121 ll'jj . the resulting configuration is then supersymmetric and has the following tension, 



where [x refers to the tension and g s is the string coupling. When we consider the physics of four 
dimensions we note that there are more general ways to get string solutions by utilizing the branes 
of string theory. If we were to wrap n-1 of the spatial dimensions of an n-brane around a compact 
cycle then such an object would appear as a string in four dimensions, with couplings and tensions 
that depended upon the details of the compactification. One model which has been worked on in 
detail is the KLMT model |14j . using a warped compactification with matter fields living on a D3 
brane in the throat. The tension of objects is then affected by the warping, with the redshift factor 
allowing for a large range of effective string tension. In principle there can be many such throats, 
each with a different warping, this would then produce strings with a variety of different tensions. 
Modelling such situations may well follow similar lines to those put forward here, adding an extra 
U(l) for each species. 

Modelling the evolution of a (p, q) network involves different approaches to that of the standard 
Nielsen-Olesen string due to the way (p, q) re-connect as two strings pass through each other. 
Abelian strings simply swap partners as they interact |15l 1161 1171 118j but this channel is not open 
to (p, q) strings as this would violate charge conservation. For such partner-swapping to take place 
an intermediate string is created, joining together the two initial strings in much the same way 
as for non- Abelian gauge vortices |19[ I20j . This raises the possibility that these networks get 
more tangled as they evolve, eventually freezing out as the Universe continues to expand. This 
would cause the string network to eventually dominate the energy density of the Universe |21| I22j . 
An important issue therefore is to decide whether these networks reach a scaling solution |23j . 
meaning that the string energy density forms a constant fraction of the background energy density, 
with the network length scale increasing proportional to time. Analytic and numerical approaches 
are now being developed to approach this problem |24l 1251 126j . here we present a simple model 
containing gauged vortices which has many of the required properties, moreover it is well suited 
to numerical implementation by adapting the available codes for evolving Abelian gauge vortices. 
This is in counter-distinction to modelling these networks with non-Abelian gauge field theories, 
where imposing lattice gauge symmetry is numerically intensive |27j . 
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The issue of the correct field theory to use when describing the low energy physics of D-strings in 
four dimensions is now converging to a supersymmetric model containing: chiral fields transforming 
in the usual way under U(l) gauge transformations; a chiral field, representing an axion, with a shift 
symmetry that gets gauged under the U(l); and a vector multiplet to carry the magnetic flux of the 
vortex[2B]. This builds on earlier work where the axion multiplet was not included |29M301 13T] . An 
obvious drawback of the model put forward here is that it is not based on a supersymmetric field 
theory. However, the simple form of coupling between the two U(l) sectors in our model leads us 
to believe that it should be possible to connect two supersymmetric U(l) gauge theories to yield a 
theory with vortex bound states of the type described below, whether such vortices would be BPS 
remains to be seen. 

We organize the paper by first introducing the model in section[2]and describing how the vacuum 
manifold depends on the parameters. The parameter space is then reduced by considering the BPS 
limit of the theory in section |2~T1 In section |3] we construct the ansatz for the vortices and calculate 
their energy, we then draw our conclusions in section^ 



2 the model 

The model is based on the simplest way to couple two independent U(l) gauge theories and has 
been seen before, albeit in a different guise jHHESl- There it was used to construct cosmic strings 
which could form a condensate in their core and become superconducting. For that to work the 
parameters are chosen such that one of the Higgs fields acquires a non-zero vacuum expectation 
value (vev) while the other sits at zero; here we require the "opposite". What we need is for the 
vacuum to be such that both Higgs fields acquire a vev, allowing for there to be independent string 
solutions following the usual arguments Although our action is the same as those for the 

superconducting string models, because we use the opposite parameter regime we write it in a way 
more suited to our analysis 

C = -D^D^-V^V^-^ U F^ -^ U T^ -VM,]^). (2.2) 

Where we have the following gauge-covariant derivatives, field strengths and potential, 

D[t = dp- ieA^, V^ = d^- igB^, (2.3) 
Fp, = dpAv-dvAp, ^ v = 8^B U - d v B^ (2.4) 

V{\<t>\M\) = y(#-77 2 ) 2 + y(^-^ 2 ) 2 -«(#-^)(^-^ 2 )- (2-5) 

For generic values of the coupling constants we find that the critical "points" are given by, 

(#,^) = (0,0), (rj 2 ,u 2 ), (r? 2 -^V,0), (Oy-^V). (2.6) 

When constructing the superconducting strings it was the last two critical points that were of 
interest, as these have one broken and one unbroken U(l) symmetry, and so the parameters were 
chosen to make these minima. As we are going to need two broken U(l)s then we need to make 
the critical point (rj 2 , v 2 ) a minimum. To gain further insight on the vacuum structure we complete 
the square as follows 

VM,m) = Mv^(#-^ 2 )±v^(#-^ 2 )l 2 -[«±iv^](#-^)W-^ 2 )(2.7) 
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and note that there are special values of k where the manifold of minima increases in dimension. 
For k = — i-s/AjA^ we get a vacuum which is an ellipsoid, 

7aT00 + y/)wH> = \AiV + V^ 2 , (2-8) 
while for k = +\\f\\\2 we get a vacuum which is a hyperboloid 

\f\[4>4> - \f\^> = \AiV - \f\~2v 1 . (2.9) 

It is these special values of k which separate the different classes of vacuum structure. For k < 
— 2\/AiA2, the minima are at (rj 2 — j^v 2 ,0), (0,z/ 2 — j^r/ 2 ), which is the parameter range for 
superconducting strings. For — ^\Ai^2 < n < +g\Ml^2j the minima are at (r] 2 ,u 2 ), which is 
the range we are interested in. For k > +^\/\i\2 the potential is unbounded below and so is 
unphysical. Within our range of interest, — ^\/Ai^2 < k < +^\Ai^2 there is another special value 
of k, namely zero. For k = we see that the two U(l) are entirely decoupled and so the vortices 
of U(l)^4 do not talk to those of XJ(1)b, in that sense we could think of two different vortex species 
as being neutrally bound. It is therefore natural to expect that k = is a boundary between the 
two vortex species having negative or positive binding energy, i.e. forming bound states or not. To 
see which sign of k gives bound states consider two separated vortices, one of each species. At the 
A vortex we have eft = 0, \ip\ J$ v while at the B vortex we tp = 0, \(f>\ i$ rj. Now consider the A 
vortex, we note from (|2.5j) that for k > we can lower the potential by reducing \tp\, this can be 
achieved by bringing the vortices closer, i.e. they attract each other. Putting together all of the 
above arguments leads us to the following parameter regime for bound A — B vortices 

< k < ^V^M- (2.10) 

2.1 BPS considerations 

Although we understand the parameter space there are still a large number of parameters, and if 
we are to make any progress we need to make some choices. To make further restrictions we take 
our inspiration from the Bogomol'nyi limit of vortices, |351 134j . Adapting those arguments to our 
case with k = the Hamiltonian takes the particular form of a sum of squares and boundary terms 
if we choose Ai = 2e 2 and A2 = 2g 2 . With that, one then sees that for a state with winding number 
n in U(l)^4 and m in U(l)# the minimum energy configuration has 

E(k = 0,n,m) = 2vrr/ 2 n + 2vri/ 2 m. (2.11) 

Where we have made the further parameter choice e = g = 1. In order to make a comparison 
between the bound-state energy of our system and (|1.1|) we introduce the following quantity, 

e M = v^vn? 2 ) 2 + (2^ 2 ) 2 . (2.12) 

What we shall now do is explore the remaining parameter space (rj, v and n) to see how the mass 
of (n,m) vortices compares to (p, q) strings. The reason we keep r/ and v is that, from (|2.11|) . they 
are expected to govern the tension of the basic constituents with \ip ~ 2irr] 2 and hf/9s ~ 2irv 2 . 

3 Vortex solutions 

What we shall now do is proceed to an evaluation of the tension of some low lying bound states. 
To be specific we shall look at two cases: firstly we consider 77 = v\ secondly we look at v = 2r] in 
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order to give a mass separation between the A and B vortices, just as fundamental and Dirichlet 
strings have different tensions. 

The ansatz we use for a static vortex with winding numbers (n, m) is as follows, 



eft = rjf(r)exp(in9), (3.13) 

Tl 

A e = -a{r), (3.14) 
e 

ip = vp(r) exp(im6>), (3.15) 

Be = ~P(r), (3.16) 
e 



for which the equations of motion following from (|2.2|) are 



f" + -/' ~ ^ f(a - l) 2 - \W{f - 1)/ + «„ V - 
r r z 2 

p" + V - ^PC9 - I) 2 - \^2v\v 2 ~ l)p + ^ 2 (/ 2 

a a — 2e f ri [o 



r 
1 



.2 2. .2/ 



-1)/ 


= o, 


(3 


-l)p 


= o, 


(3 


-1) 


= o, 


(3 


- 1) 


= 0. 


(3 



Note that we have kept e, g, \%, A2 in these equations even though we shall only be using the 
values as described above, e = g = = 5 A2 = 1. The asymptotic boundary conditions are that 
f(r — * 00) —* 1, p(r — > 00) — > 1, a(r — > 00) — > 1, /3(r — > cxd) — » 1, necessary for a finite energy 
configuration. The boundary conditions at the vortex core depend on the winding number; if <fr 
has a non-zero winding then /(0) = and a(0) = otherwise f'(0) = and a(0) = 1. Similar 
conditions hold for ip. 



3.1 Case I, 77 = v 

The first example is for rj = v and we can use our understanding of the potential gained earlier 
to see how the profile functions should behave. In Fig. ^]we have a representation of the vacuum 
with: the points (±1, ±1) and (±1, =pl) being the vacua for k in our region of interest, —1 < k < 1; 
the circle is the vacuum for k = — 1; and the 45° lines are the vacuum for k = 1. On the diagram 
we have shown how the Higgs field profile functions vary as we move from one vacuum at x — » — 00 
to another at x — > 00 for a vortex with winding in the <j) field but not ip. The various solid lines 
correspond to the vortex solutions for different values of k, the straight line at ip/v = 1 is the 
k = case where <p and ip are decoupled. The lines below this case are for k increasing, and they 
approach the 45° lines as k nears unity. If we decrease k below zero then the Higgs profile functions 
mark out a curve in <f> — ip space which tends toward the circle, which is the vacuum for k = — 1. 

In order to solve (|3. 17H3.2T7^ we used a numerical approach described in jHS] called successive 
over-relaxation. We tested our numerics against the analytic BPS solution and we are confident 
that the energies quoted in the appendices are correct to four decimal places. What is clear from 
table 1 is that vortices with higher winding numbers do indeed form bound states, with the bound 
state energies being lower than the sum of their constituents. Moreover, as k is increased toward 
unity (the upper bound on k) the binding of vortices is stronger. The results of table 1 have been 
plotted in Fig. [21 with the energies scaled in order to make comparison with ()2.12j) . which has been 
shown with dotted lines. The plot consists of three "rays" corresponding to different values of m, 
with the lowest such ray being m = and the upper being m = 2. Within each ray we show 
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Figure 1: A representation of the field space, with the vortex profile functions for a (1,0) vortex 
marked. The circle is the k = — 1 vacuum and the 45° lines are the k = 1 vacuum. 

how the tension depends on k, with the upper line in each ray being k = 0.1 and the lowest being 
k = 0.9. We can see from this plot that the bound state spectrum does not follow the (p, q) tension 
spectrum of 1)2. 12j) particularly well, but for k near its upper bound it is closest. In a real network 
we would expect there to be fewer high winding strings as they are heavier, as such, taking k close 
to unity may well give a suitable approximation. 

3.2 Case II, 2rj = u 

In order to represent the difference in tension between a D-string and and F-string we can choose to 
give different values to r\ and v. For the purposes of illustration we took v 2 = Arj 2 and repeated the 
above process to produce table 2, which contains the tension of the vortices for the first few windings. 
Again we see that higher winding strings are bound states, with an energy lower than the sum of 
their constituent parts. Fig. |3] illustrates the data in table 2, with the dotted lines representing 
()2.12)) . We note that the specific tension formula for (p, q) strings seems to be modelled better in 
this second case, where the basic strings have different tensions, for k close to its upper bound. 
This means that simulations using different values for rj and v will more closely model the dynamics 
of a (p, q) network. 

4 Conclusion 

In trying to model the dynamics of a network of (p, q) strings coming from string theory one must 
try both analytic and numerical approaches (2HES1I2S]) each with its advantages and disadvantages. 
The properties of inter-commutation of (p, q) strings, i.e. the presence of three-string vertices, makes 
it natural model these networks with a non-Abelian gauge theory containing a rich spectrum of 
vortices. Unfortunately such simulations are extremely demanding in terms of cpu hours, making 
large scale simulations intractable. As one needs a large dynamic range in order to observe scaling 
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12 3 4 

n 

Figure 2: string tension for the case rf = v 2 , k=0.1, 0.2, 0.5, 0.9. Of the three "rays", the lower 
corresponds to m = 0, the middle has m = 1 and the upper has m = 2. E re f is a reference energy, 
E ref (m = 0) = £(l, ), £re/0™ > 0) = E( 0j l). 




12 3 4 



n 

Figure 3: string tension for the case 4r/ 2 = z^ 2 , k=0.1, 0.2, 0.5, 0.9. Of the three "rays", the lower 
corresponds to m = 0, the middle has m = 1 and the upper has m = 2. E re f is a reference energy, 
E ref (m = 0) = £(i,o), E re f(m > 0) = £( ,i)- 
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behaviour this presents us with a problem; what we have presented here is a practical resolution 
to this impasse. By using two Abelian gauge symmetries it is a simple matter to adapt existing 
cosmic string code to the study of (p, q) networks. 

An important matter that this work does not reveal is the details of the interaction between 
two vortices of different type. In particular we have not shown that such vortices inter-commute 
to create string junctions; it is possible that they simply pass through each other as they would for 
k = 0. However, as string junctions certainly exist in this model we expect inter-commutation to 
occur. 

A natural follow-up of this work would be to try to extend the proposed model incorporating 
the ideas of [2S]. In j^H] it was argued that D-strings can be modelled by a U(l) super symmetric 
gauge field theory containing a chiral axion field as well as the usual chiral Higgs fields. If we are 
to have a field theory description of (p, q) strings along the lines presented here one would need 
to couple together two such field theories. It is certainly possible to couple two such field theories 
together, however it is unclear that one would find supersymmetric bound-states of vortices as a 
solution. 
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K 


n 


m 


^(n,m) 


nE, 10) + m&oi) 


0.1 


1 





0.9985 


0.9985 


0.1 


2 





1.9956 


1.9952 


0.1 


3 





2.9922 


2.9928 


0.1 


4 





3.9886 


3.9904 


0.1 


1 


1 


1.9568 


1.9952 


0.1 


2 


1 


2.9321 


2.9928 


0.1 


3 


1 


3.9160 


3.9904 


0.1 


2 


2 


3.8875 


3.9904 


0.2 


1 





0.9939 


0.9939 


0.2 


2 





1.9824 


1.9878 


0.2 


3 





2.9688 


2.9817 


0.2 


4 





3.9540 


3.9756 


0.2 


1 


1 


1.9098 


1.9878 


0.2 


2 


1 


2.8565 


2.9817 


0.2 


3 


1 


3.8184 


3.9756 


0.2 


2 


2 


3.7661 


3.9756 


0.5 


1 





0.9586 


0.9586 


0.5 


2 





1.8816 


1.9172 


0.5 


3 





2.7910 


2.8758 


0.5 


4 





3.6931 


3.8344 


0.5 


1 


1 


1.7358 


1.9172 


0.5 


2 


1 


2.5677 


2.8758 


0.5 


3 


1 


3.4233 


3.8344 


0.5 


2 


2 


3.3243 


3.8344 


0.9 


1 





0.7931 


0.7931 


0.9 


2 





1.4354 


1.5862 


0.9 


3 





2.0274 


2.3793 


0.9 


4 





2.5928 


3.1724 


0.9 


1 


1 


1.2780 


1.5862 


0.9 


2 


1 


1.7981 


2.3793 


0.9 


3 


1 


2.3235 


3.1724 


0.9 


2 


2 


2.2284 


3.1724 



Table 1: Table showing the energy of the vortices for v = rj. 
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K 


n 


m 


Z71 


nE (1 o) + mEiQ ±) 


0.1 


1 





0.9981 


0.9981 


0.1 


2 





1.9949 


1.9962 


0.1 


3 





2.9912 


2.9943 


0.1 


4 





3.9874 


3.9924 


0.1 





1 


3.9960 


3.9960 


0.1 


1 


1 


4.9275 


4.9941 


0.1 


2 


1 


5.9049 


5.9922 


0.1 


3 


1 


6.8949 


6.9903 


0.1 





2 


7.9887 


7.992 


0.1 


1 


2 


8.8697 


8.9901 


0.1 


2 


2 


9.8226 


9.9882 


0.2 


1 





0.9924 


0.9924 


0.2 


2 





1.9794 


1.9848 


0.2 


3 





2.9646 


2.9772 


0.2 


4 





3.9490 


3.9696 


0.2 





1 


3.9845 


3.9845 


0.2 


1 


1 


4.8498 


4.9769 


0.2 


2 


1 


5.7993 


5.9693 


0.2 


3 


1 


6.7722 


6.9617 


0.2 





2 


7.9524 


7.969 


0.2 


1 


2 


8.7291 


8.9614 


0.2 


2 


2 


9.6335 


9.9538 


0.5 


1 





0.9480 


0.9480 


0.5 


2 





1.8598 


1.896 


0.5 


3 





2.7605 


2.844 


0.5 


4 





3.6556 


3.792 


0.5 





1 


3.9022 


3.9022 


0.5 


1 


1 


4.5727 


4.8682 


0.5 


2 


1 


5.4033 


5.7718 


0.5 


3 


1 


6.2771 


6.7462 


0.5 





2 


7.7105 


7.8044 


0.5 


1 


2 


8.2279 


8.7524 


0.5 


2 


2 


8.9735 


9.7004 


0.9 


1 





0.7388 


0.7388 


0.9 


2 





1.3294 


1.4776 


0.9 


3 





1.8811 


2.2164 


0.9 


4 





2.4134 


2.9552 


0.9 





1 


3.5994 


3.5994 


0.9 


1 


1 


3.9531 


4.3382 


0.9 


2 


1 


4.4594 


5.077 


0.9 


3 


1 


4.9896 


5.8158 


0.9 





2 


6.9420 


7.988 


0.9 


1 


2 


7.1733 


7.9376 


0.9 


2 


2 


7.6053 


9.4656 



Table 2: Table showing the energy of the vortices for v = 2rj. 
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